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Abstract 

There exists a large number of experimental and theoretical results supporting the picture of "macroscopic 
qubits" implemented by nanoscopic Josephson junctions. On the other hand the standard model of such systems 
given in terms of a single degree of freedom suggests their semiclassical behavior due to a localization mechanism 
caused by a strong coupling to an environment. Indeed, such a mechanism is observed in an atomic Bose-Einstein 
condensate (BEC) placed in a double- well potential - a system mathematically equivalent to a Josephson junction. 
In this note it is shown, on the example of a Cooper pair box, that replacing the BEC-type model for Cooper 
pairs by a lattice gas model one can reduce the environmental effects of "dequantization" and can explain the 
experimental data in particular the existing huge differences between the measured values of relaxation times. 

1 Introduction 

In the last decade remarkable experiments were performed involving measurements and manipulations of states for 
a single or several nanoscopic Josephson junctions which were consistently interpreted in terms of two level quantum 
systems [TJ El El [U \E\ . The main assumption in the theoretical analysis is that such a many-body mesoscopic system 
can be effectively treated as a quantum system of a single degree of freedom typically described by a large spin 
or nonlinear oscillator model. A nonlinear Hamiltonian yields the structure of two lowest energy levels which at 
the low enough temperatures can be separated from the others to form an effective macroscopic qubit. The main 
problem with such models is a presence of a typically strong and collective coupling to an environment. Namely, it 
is expected that the observed states should be rather well-localized semiclassical ones, which seem to be the only 
relatively stable with respect to external noise. This mechanism in briefly discussed in the next Section. However, 
the semiclassical states for the model of Cooper pair box (CPB) are characterized by large charge fluctuations 
which are not observed in the experiments. Therefore, either environmental decoherence producing semiclassical 
states does not work for Josephson qubits at the typical time scale of the experiments or the standard single degree 
of freedom model is not correct. The first alternative seems to be unlikely because the semiclassical character of 
observed states is confirmed in the recent experiments on atomic Bose-Einstein condensate (BEC) in a double-well 
potential [7]. Although this is a physically different system its mathematical description is the same as for the 
standard model of a CPB. In the Section 3 the second alternative is discussed, a lattice gas model of CPB. In this 
model the coupling to an environment has individual and local character and the effective size of the system is much 
smaller than in the standard model. Therefore, the picture of an effective qubit can be correct and some quantum 
features can be present on the relevant time scale. The price we pay is that the qubit states are not as well defined 
as in the standard model and strong leakage processes to other states in the effective Hilbert space are present. 



2 Collective single degree of freedom model 

The most studied, both experimentally and theoretically, examples of mesoscopic system which should support a 
qubit are "superconducting qubits". For simplicity, only a CPB called "charge qubit" is discussed here. It is a 
circuit consisting of a small superconducting island connected via Josephson junction to a large superconducting 
reservoir. Coulomb repulsion between Cooper pairs in a small electrode become important and must be taken into 
account in the Hamiltonian. The simple Josephson Hamiltonian reads pQ 

H = AE C £(n - n - n g ) 2 \n)(n\ - ^ £(|n + l)(n| + \n)(n + 1|) (1) 
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where |n) describe the state with n Cooper pairs on the island, Eq determines the magnitude of the Coulomb 
repulsion, Ej governs the tunneling process, no >> 1 is a number of Cooper pairs on the island at the neutral 
reference state and the additional n g is a fine tuning of the external control. To establish a typical value of no is 
a subtle problem, some authors assume that all Cooper pairs should be taken into account what yields n ~ 10 8 , 
on the other hand one can argue that only the electrons close to Fermi surface should matter what gives n ~ 10 4 . 
The main assumption behind the simple Hamiltonian {TJ is that all Cooper pairs occupy a single quantum state, 
similarly to the BEC picture and therefore only the number of Cooper pairs matters. 

Under the assumption no >> 1 and restricting to the states with |n — no| << no the Hamiltonian ([I]) can be 
rewritten in terms of spin variables Jk, k = 1, 2, 3 satisfying standard relations 
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[Jk, Jl]=iJ2 e klmJm , Jl + J% + Jl = j(j + 1) (2) 
m— 1 

with j = no. The Josephson Hamiltonian reads now 

H = 4E c (J 3 -n g ) 2 -f^Ji. (3) 
An essentially equivalent model Hamiltonian has a form of the nonlinear oscillator's one 

H = 4E c (ala-(n +n g )) 2 --^(a + al) (4) 



where [a,a^] — 1. In the following only the large spin model is used. 

The device is controlled by external electromagnetic fields which are coupled to the net electric charge Q = 2e J3 
and to the electric current = i[H, Q] ~ J2. Hence the control Hamiltonian is given by 

H c (t) = h 2 (t)J 2 + h 3 (t)J 3 . (5) 

To describe the influence of an environment one should notice that the leading contribution to the system- 
bath interaction is always of the form similar to (J5j) with the external fields hk(t) replaced by bath's operators 
Bk] k — 1,2,3 and taking into account that the fluctuations of the tunneling rate produce the term with k = 1. 

Under quite general conditions the reduced dynamics of the spin in Markovian approximation [8] is given by the 
master equation for the reduced density matrix p 



d 

-p = -i[H,p} + J2 c k i([JkpJi} + [Jk,pJi}) (6) 

k,l=l 



where H is given by ([3]) and the positively defined matrix [cm] depends on the details of the reservoir. The stability 
of the initial pure state can be characterized by the initial decay of purity 

d 3 

r(V0 = --TrpVo = 2 Y, C «((V^ JJnP) ~ (V^KM^)) (7) 

k,l = l 

The RHS of (7|) can be treated as a dispersion of the total spin modified by some weights 7^- eigenvalues of 
the relaxation matrix [cm]- Therefore, one can see that the purity of the eigenstates \j,m) decays with the rate 
r(|j, to)) ~ "fj 2 . By semiclassical states we mean the states which minimize the dispersion of the total spin to the 
value of the order ~ j. The decay rate of such states is of the order of *yj (7 is a typical value of the relaxation 
constants 7^). It follows that for large enough spins only semiclassical states are relatively stable and can be 
observed [6]. 

This type of "dequantization" is supposed to be the main mechanism of emergence of classical world and 
provides the solution for the problem of molecular structure, existence of deformed nuclei and generally the absence 
of Schrodinger cat states. A nice illustration of this mechanism is provided by the experiments with atomic BEC 
placed in a double- well potential • The Hamiltonian of such a system can be approximated by the Josephson one 
([3]) where the observable J3 corresponds to the excess number of atoms with respect to an equilibrium value no in 
a chosen well. The semiclassical fluctuations of the number of bosons in a given well ~ ^/no are clearly observed. 

On the other hand, for CPBs such large fluctuations of the number of Cooper pairs are not observed, on the 
contrary, superpositions of states which differ by a single Cooper pair are visible in numerous experiments. This 
suggests that perhaps the model based on the Josephson Hamiltonian |(3]) is not adequate. 
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3 Lattice gas model 



In this model one assumes that due to strong Coulomb repulsion between Cooper pairs on the small electrode they 
cannot occupy a single quantum state but rather a set of localized sites with the occupation number 1 or 0. The 
tunneling to a large electrode is possible only from the sites placed at the boundary of the small electrode which 
is close to the large one. The number of these relevant sites, denoted by m should be roughly proportional to the 
square root of the total number of Cooper pairs n . Taking n ~ 10 8 — 10 4 one obtains m ~ 10 4 — 10 2 . Introducing 
the occupation number operators hk — + l)/2; k = 1,2, mo with eigenvalues 0, 1 and spin-1/2 Pauli matrices 
<7fc, a — x,y, z, +, — , k = 1, 2, mo one can propose the following model Hamiltonian for this system 

too p mo 

H = 4£ C QT n k - (m + n g )f - -± + (8) 

fc=l k=l 

Here the first term describes Coulomb repulsion and n g 6 [0, 1) is a fine tuning control parameter. The second term 
describes tunneling from and to individual sites with the normalization of local tunneling parameters X) ICfc| 2 = 1. 
Using the standard assumption Ec > Ej » kT one can restrict the effective Hilbert space to the lowest energy 
spectrum sector 7i e s spanned by the following m + 1 basis eigenvectors of the operators fik 

|0) = , \k) = |1,...,1,0,1,...,1) , fc=l,2,...,m . (9) 

k 

Introducing a normalized vector and two projectors in 7i c a 

mo 

I0 = 5>l*> - p = lo)(o| + |0(£l , P ± = i-P (io) 
fc=i 

one can write down the Hamiltonian ([8]) restricted to H c s and up to an irrelevant constant as 



H = ^(IO<£l " I0>(0| - P) - ^(|0(0| - I"-:.; ' ■ = 4£< (1-2,,,,). (11) 



The energy levels of ifTTj) consist of the mo — 1 manifold of degenerated levels of the energy equal to zero given by 
the projector P and two levels |±) of positive and negative energies E± 

|+>=cos~|0-sin^|0) , |-) = cos ||0) + sin 6 -\^) , E± = ±\<J E(n B ) 2 + Ej - ~E(n g ) (12) 



^ = |+X-| ,&7 = (a+)t , af = (|+X+| - |-)(-|). (14) 



where 9 is defined by cos 9 = E(n g ) j ' yj 'E(n g ) 2 + E'j. 

The external control is performed by the coupling through the total electric charge operator Q and the total 
electric current J which, when restricted to 7i g reads 

Q = e(|0(£|-|0)<0|+P x ) , J = i[H,Q}=ieEj(\Q){^-\0{0\). (13) 

Obviously, if the system is completely isolated the qubit space spanned by |0) and |^) is invariant with respect to 
the Hamiltonian and the external control yielding the usual model of charge qubit. We denote its qubit observables 
by 

In contrast to the large spin model the coupling of the lattice gas model to a bath has individual and local 
character. The sufficiently general form of the system-bath interaction Hamiltonian is given by 

toq 

Hint=J2 E **®£fc ( 15 ) 
fc=l a=x,y,z 

restricted to the subspace H c s ® Hbath- One can assume that the bath operators B% correspond to different 
independent constituents of the bath ("private baths" picture). Applying now the standard weak coupling limit 
technique and assuming that the temperature of the bath is zero and the private baths are identical one obtains 
the following master equation for the reduced density matrix of the system [8] 

j t P - - i[H, p] + \([a-p, ct+] + [a-,pa+]) ~ z#f , [«r|, p\\ (16) 

mo — 1 mo — 1 

+ PiJ2 {[a k p,al} + [a k ,pal])+p j2 £ ([hp, b{] + [b k , pb\}). (17) 

k=l k=l 
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Here 

& k = P ± \k)(+\ , b k = \-)(k\P x (18) 

and the decay rates A, /Ui,/i2 > and "pure dephasing" rate v are given by standard expressions involving Fourier 
transforms of the autocorrelation functions for baths operators taken at frequencies corresponding to proper energy 
differences. 

One can obtain from (fT?]) the Pauli Master Equation for the level occupation probabilities p± — (±|/5|±) , 
p = Tr(pP^): 

^ = -2[A+(m -l)/xi]p+ (19) 
— = 2(m - l)pip+ - 2^ 2 Po (20) 
^ = 2Ap+ + 2/i 2 p (21) 
and the closed equation for the qubit coherence z = Tr(paf ) 



dz I ' 

— = (iu - [A + Av + (m - l)m])z , oj = E{n g ) 2 + E 2 , (22) 



One can notice that the leading damping effects in l|2H |22|) are proportional to T = (mo — l)pi >> [1%, /i 2 , v, A. 
Taking only them into account one obtains an approximative solution on the time scale of the order T _1 

p+(t) = P+ (0)e- rt , p_(t)=p_(0) , po(t) =Po(0)+p+(0)[l-e- rt ] , z(t) = z(0) exp{(z^ - T/2)t}. (23) 

In this approximation the third, highly degenerated level acts a "probability sink" for the qubit. The long time 
behavior is given by 

, P+ (t) = 0, p (t)=p+(0)e- 2 ^ t , p_(t) = l-po(i) , «(t) = 0. (24) 

The only directly measured observable is the net charge 

Q(t) = Tr(p(t)Q) = (z{t) +z(jt))em6- (1 - p (t)) cos 6 +p (t). (25) 

We are interested also in the time dependence of the averaged energy given by 

£(t) = Tr(p(t)H) = E +P+ {t) + B_p_(t). (26) 

Assume, for simplicity, that one prepares a system in a pure state |£) and choose n g — 1/2 (sm6 = 1) what 
corresponds to the experimental setting of [2], j3|. Combining the formulas 1)251) and f24|) one obtains the evolution 
of the mean charge for short times 
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Q(t)= (cos(urf) + -e-~ e-~ + -. (27) 
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Although the predicted form of "coherent charge oscillations" differs from that obtained for the 2-level model 

Q(t) = (2cosLut)e~T5 , (28) 

with the decoherence time T 2 , both curves can be fitted to the experimental data. Moreover, both curves give the 
same Lorentzian frequency spectrum around u when T 2 = 2/T. 

The measurements of energy relaxation are usually performed for n g ~ (sin 9 ~ 0) . As Ec > Ej hence 

E+ ~ , E^ ~ -4E C ■ (29) 

It follows from l(2~4|) . l|26p and (|29|) that the energy decays to its lowest value on the long time scale 

£(t) =4E c (p+(0)e- 2 ^ t -l). (30) 

The relation (|30|) can explain observed in [4] long relaxation times corresponding to Ti = (2^ 2 ) _1 in contrast to 
a much faster decay of p+(t) with decay time T\ = T^ 1 which seems to determine the results of [3]. The ratio 
T 1 /T 2 ~ m /2 estimated from the experimental data of [4] gives a reasonable value m — 10 3 . 
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4 Conclusions 



A new model of a Cooper pair box is proposed, based on the picture of a large number too of localized states 
("lattice sites") which are occupied by at most one Cooper pair for each site and from which the Cooper pairs can 
tunnel to a large superconducting reservoir and back. This model is compared with the standard model of a large 
spin j = n system, where n is an average number of Cooper pairs which occupy a single quantum state delocalized 
over the whole superconducting island. The main difference between both models is due to the interaction with an 
environment. For the large spin model all Cooper pairs interact collectively with the bath what makes the states 
with a fixed charge very unstable with decoherence rates <~ rig. The relatively stable states with life-times <~ tiq 
possess semiclassical character and yield large charge fluctuations of the order ^Juq. This behavior is confirmed for 
the atomic Bose-Einstein condensate in a double well, which is described by the similar mathematical model. It 
is very difficult to explain why this mechanism should not work for a CPB and why the well defined charge states 
remains relatively stable. 

The lattice gas model suggests a localized coupling of any site to its individual heat bath what at the low 
temperature regime produces charge fluctuations of the order 0(1) in agreement with experiments. However, the 
many-body character of the system is still present. The relevant part of the Hilbert space is spanned not only 
by a pair of Hamiltonian eigenstates corresponding to a qubit but also by to — 1 degenerated states which form 
a probability sink for the excited state of the qubit with the decay rate ~ (to — 1). The predicted form of the 
damped charge oscillations differs from the standard formula providing a possible experimental test of the model. 
This model explains also a huge discrepancy between different data on relaxation times, introducing natural two 
time scales for dissipative processes in the system. 

Finally, one should mention another difference between the standard qubit model of CPB and the lattice gas 
model. In the former the effective 2-dimensional Hilbert space is spanned by the vectors corresponding to non- 
degenerated eigenstates of a charge operator. In the later model, one of the qubit states is a superposition of 
a large number of degenerated charge eigenstates with the probability amplitudes determined by local tunneling 
rates which can vary in time. Therefore, one can expect to observe a slow random drift of that state over the mo 
dimensional Hilbert subspace. This means that not only the probability leaks from the qubit but the qubit itself 
does not correspond to a fixed well-defined 2-dimensional subspace. 
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